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Abstract 

The earthworm family Hormogastridae shows a remarkable disjunction in its distribution in the Iberian 
Peninsula, with the Hormogaster elisae species complex isolated from the rest of the species. Hormogaster 
joseantonioi sp. n., a new species found in the intermediate area berween the main ranges (in Teruel, 
Aragon), was described following the integrative approach, as it is suitable for earthworms due to their 
highly homoplasic morphology. The phylogenetic analysis of the molecular markers placed the new spe- 
cies as a sister taxon to H. elisae, thus showing the colonizing lineage of Central Iberian Peninsula could 
have originated near the H. joseantonioi sp. n. current range. External motphological characters revealed 
some degree of overlap with previously described species, but internal characters presented configurations/ 
states unknown from other members of the family. These traits make the new species a key piece to un- 
derstand the evolution of Hormogastridae. 



Copyright D.E Marchan et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), 
which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 
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Introduction 

The increasing availability of molecular and ecological data has placed the integrative 
taxonomy (as defined by Dayrat 2005) as a viable alternative to traditional species 
description. Several authors advocate its use in different animal groups (Padial and De 
La Pviba 2010; Schlick-Steiner et al. 2010; Heethoff et al. 201 1; but see Yeates et al. 
2011 for iterative taxonomy instead) and particularly in earthworms (Blakemore and 
Kupriyanova 2010; Novo et al. 2012), whose taxonomy is in need of deep revision in 
the light of molecular phylogeny (Jamieson et al. 2002; Pop et al. 2003, 2007; Chang 
et al. 2008; Briones et al. 2009; Perez-Losada et al. 2009, 2011; Novo et al. 2011; 
Fernandez et al. 2012). 

Fernandez et al. (2014) have developed a new tool based in micro-computed to- 
mography to study specimens in a non-destructive way which could help as an addi- 
tional source of information. 

Taxonomic characters traditionally used for the study of earthworms are few and 
sometimes present high intraspecific variability (Michaelsen 1900 and Stephenson 
1930 on their global fauna; Pop et al. 2003 and Briones et al. 2009 about lumbricid 
earthworms). Recent findings show that cryptic diversity is common in these animals 
(but see critique in Blakemore et al. 2010), therefore earthworm taxonomy can par- 
ticularly benefit from an integrative approach. 

Novo et al. (2011) presented a molecular phylogeny of Hormogastridae (Oli- 
gochaeta, Annelida), whose taxonomy has historically been built on morphological 
characters, which highlighted some interesting evolutionary aspects. On one hand, 
hormogastrid distribution across the Western Mediterranean is biogeographically con- 
sistent, reflecting the geological events that affected the region in the Tertiary (which 
confirms previous studies, e.g Bouche 1972, Sbordoni et al. 1992). Two species -Xana 
omodeoi Diaz Cosin, Briones & Trigo, 1989 and the morphospecies Hormogaster eli- 
sae Alvarez, 1977 -, however, are found in locations far apart from the family main 
range in the Iberian Peninsula. While all the other Iberian species are distributed in 
Northeastern Spain, X. omodeoi inhabits Northwestern Spain and H. elisae is found in 
Central Spain (Segovia, Madrid and Guadalajara). The result is a disjunct distribution. 

Novo (20 1 0) found Hormogaster elisae complex to be monophyletic, and thus the 
likely result of a single colonisation event presumably from the North or the East of 
the Iberian Peninsula. There could be remaining populations of the migrating lineage 
in the geographic gap, which haven't been discovered yet. 

On the other hand, it seems that most key characters used for hormogastrid tra- 
ditional taxonomy and phylogeny (notably the shape, number and position of the 
spermathecae) are highly homoplasic, showing little or no phylogenetic signal across 
the family. 
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Due to its relevance for this subject, the intermediate area between the main rang- 
es of hormogastrids in Spain has been subject to recent sampling campaigns. Both 
Zaragoza and Teruel (Aragon, Spain) were suitable regions as they have been poorly 
sampled for earthworms unlike the surrounding provinces. While no success was met 
in Zaragoza, a population assignable to a new species of Hormogastridae was recently 
found in Teruel. 

This paper focuses on the description of Hormogaster joseantonioi sp. n. from an in- 
tegrative point of view, following the example of Novo et al. (2012). The new molecular 
and morphological data are interpreted to gain insight into the diversification and mor- 
phological radiation of the family, with some considerations about its constituent genera. 



Materials and methods 

Earthworm specimens and sampling points 

Specimens were collected by hand and fixed in the field in ca. 96% EtOH, with sub- 
sequent alcohol changes. Once in the laboratory, specimens were preserved at -20 °C. 

The studied material includes 1 0 specimens (five mature specimens, one semima- 
ture specimen with tubercula pubertatis and four immatures) collected in a cleared 
holm-oak wood at the foothill of Sierra de Oriche, road A-2514 between Huesa del 
Comun and Rudilla, Teruel (Spain) (4l°0'55.68"N, 0°58'55.98"W) (Figure 1). 




Figure I . Map of the Iberian Peninsula showing the collection site of H. joseantonioi sp.n. (indicated by 
the white star). The northeastern hormogastrid range is shown in green, H. elisae range is shown in pink 
and X. omodeoi known location is indicated in yellow. 
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Specimens have been deposited in the Oligochaete collection of the Departamento 
de Zoologia y Antropologia Fi'sica, Universidad Complutense de Madrid (UCMLT), 
Spain with vouchers UCMLT 00001-00010. 

Specimens available from previous studies (Novo et al. 2010, 2011, 2012) of all 
known hormogastrid species were used for comparison. Morphological characters in- 
clude those features traditionally used for hormogastrids and other earthworms. 



Molecular data generation 

Total genomic DNA was extracted from ventral integument tissue samples using the 
DNeasy Tissue Kit (QIAGEN) with two consecutive steps of elution (70 ul of buffer). 
Seven molecular regions were amplified: mitochondrial subunit 1 of cytochrome c oxidase 
(COI), 16S rRNA and tRNA Leu, Ala, and Ser (16S t-RNAs), one nuclear ribosomal 
gene (a fragment of 28S rRNA) and one nuclear protein-encoding gene (histone H3). 
Primer sequences, polymerase chain reactions (PCR) and sequencing reactions are the 
same as in Novo et al. (201 1). GeneBank accession numbers for the holo- and paragene- 
types, following Chakrabarty (2010) for the markers analysed here are shown in Table 1. 



Phylogenetic analyses 

The new sequences were combined with all the hormogastrid information available 
from previous studies (Novo et al. 2010, 2011, 2012) in order to find their phylo- 
genetic placement inside the family Pontodrilus litoralis Grube, 1855, Dichogaster 
saliens Beddard, 1893, Amynthas robustus Perrier, 1872, Lumbricus terrestris Linnaeus, 
1758 and Aporrectodea trapezoides Duges, 1828 were used as outgroups (all the Gen- 
Bank accession numbers are shown in Appendix). As hormogastrid individuals from 



Table I . Holo- and paragenetypes (sensu Chakrabarty, 2010) of H. joseantonioi sp. n., and their Gen- 
Bank accession numbers. The hologenetype is shown in bold. 



Specimen 


Voucher 


COI 


16S-tRNAs 


28S rRNA 


H3 


HRUD1 


UCMLT 00001 


KJ632674 


KJ632684 


KJ632686 


KJ632688 


HRUD2 


UCMLT 00002 


KJ632675 


KJ632685 


KJ632687 


KJ632689 


HRUD3 


UCMLT 00003 


KJ632676 








HRUD4 


UCMLT 00004 


KJ632677 








HRUD5 


UCMLT 00005 


KJ632678 








HRUD6 


UCMLT 00006 


KJ632679 








HRUD7 


UCMLT 00007 


KJ632680 








HRUD8 


UCMLT 00008 


KJ632681 








HRUD9 


UCMLT 00009 


KJ632682 








HRUD10 


UCMLT 00010 


KJ632683 
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the same locality usually cluster together, one individual was analysed as representa- 
tive per sampling site. 

Sequences of each individual gene were aligned in MAFFT (Katoh and Standley 
2013) with default settings and concatenated, resulting in a matrix of 2532 bp. jMod- 
elTest v. 2.1.3 (Darriba et al. 2012) was used to select the best-fit evolutionary model 
using the Akaike information criterion (AIC; Akaike 1973), and Bayesian informa- 
tion criterion (BIC; Schwarz 1978) which were GTR+I+G for COI, 16s and 28s, and 
HKY+I+G for H3. 

Bayesian Inference (BI) of the phylogeny was estimated with MRBAYES v.3.1.2 
(Ronquist and Huelsenbeck 2003) implemented in the CIPRES Science Gateway 
V. 3.3. (http://www.phylo.org/index.php/portal/). Unlinked nucleotide substitution 
models selected were specified for each gene fragment and the nucleotide substitu- 
tion estimates were allowed to vary independently between each partition. Param- 
eters were set to ten million generations and 10,000 trees were sampled for every 
1000th generation, initiating the analysis from a random tree. After two analysis 
were performed 20% of the trees were discarded as burn-in, and the remaining trees 
were combined to find the maximum a posteriori probability estimate of phylogeny. 
Maximum likelihood analyses were performed with RAxML 7.2.7 (Stamatakis 2006) 
in the CIPRES Science Gateway with default settings, using GTR+I+G for each data 
partition and estimating the support for the resulting topologies by 100 bootstrap 
replicates. 

Uncorrected pairwise differences for the mitochondrial regions were calculated be- 
tween H. joseantonioi and the most closely related species with Arlequin 3.5 (Excoffier 
and Lischer 2010. To visualize the genetic distance we constructed networks with 
SplitsTree4 v.4.1 1.3 (Huson and Bryant 2006) for the more variable genes, including 
the former species plus Hormogaster riojana Qiu & Bouche, 1998 and A. trapezoides as 
outgroups. Default settings were used. 



Results 

Taxonomic results 

Phylum Annelida Lamarck, 1802 
Subphylum Clitellata Michaelsen, 1919 
Class Oligochaeta Grube, 1850 
Superorder Megadrilacea Benham, 1890 
Order Haplotaxida Michaelsen, 1900 
Family Hormogastridae Michaelsen, 1900 

Genus Hormogaster Rosa, 1887 
Type-species. Hormogaster redii Rosa, 1887. 
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Hormogaster joseantonioi Fernandez Marchan, sp. n. 

http://zoobank.org/lB7B13C0-FA56-466E-9FFE-AB985EB582BA 
http://species-id.net/wiki/Hormogaster_joseantonioi 

Material examined.Holotype.Adult(UCMLT00003),4l°0'55.68"N,0°58'55.98"W, 
from a cleared holm-oak wood on the foothill of Oriche mountains, road A-2514 be- 
tween Huesa del Comun and Rudilla, Teruel (Spain), collectors D. Fernandez Marchan 
and JA. Fernandez Fernandez. 

Paratypes. Nine individuals (UCMLT 00001, 00002, 00004-00010), with the 
same collection data of the holotype. 

Other material examined. 1 6 hormogastrid species and several subspecies belong- 
ing to the UCMLT collection. 

Morphological description. External morphology (Figure 2). *Measures taken on 
the two only complete specimens, one being the holotype. 

Length of mature specimens*: 178—180 mm. 

Maximum diameter (pre-clitellar, clitellar, post-clitellar) of mature specimens: 
8-10, 9-11, 7-10 mm. 

Number of segments*: 305-369. 
Weight (fixed specimens)*: 7.05-1 1.57 g. 

Colour: From light brown to dark chocolate brown varying between individuals, 
with orangeish-brown clitellum of a lighter shade on living specimens (Figure 2a). 
Beige with brown stripes or patches, mainly on the anterior end, with darker clitellum 
on fixed specimens (Figure 2b). 

Prostomium prolobic, longitudinal striation on segments 1 and 2. 

Closely paired chaetae; interchaetal ratio at segment 40, aa: 33, ab: 1.3, be: 6, cd: 
1, dd: 27. Nephridial pores in a row between chaetae b and c (very close to b), visible 
on fixed specimens as a brownish line. 

Spermathecal pores at intersegments 9/10 and 10/11 at the level of cd. 

Male pores open over chaetae ab at the intersegment 15/16, surrounded by heart- 
shaped porophores which extend over most of segment 15 and at least half of 16. 
Female pores in segment 14 at the same level as male pores. 

Clitellum saddle-shaped extending over segments (13) 14-28. Tubercula puberta- 
tis on 1/n 22-27(1 /n 28) as a continuous line. Papillae of chaetae ab in variable posi- 
tions, usually between segments 12 and 28: papillae on 12 always showing an unusual 
degree of development in mature individuals, being very conspicuous both in live and 
fixed specimens (Figure 2a). 

Internal anatomy. Funnel shaped, strongly thickened septa in 6/7, 7/8 and 8/9, 
septum 9/10 slightly thickened. The latter's attachment to the dorsal body wall is 
displaced two segments backwards, creating a mismatch between inner and outer seg- 
mentation with an internally very wide segment 9. 

Last pair of hearts in segment 1 1. Three shiny, strongly muscular gizzards in 6, 7 
and 8. Not apparent Morren's glands, even though small wrinkles exist in the oesopha- 
geal wall between segments 10 and 16. 
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Figure 2. (A) Live specimens of H. joseantonioi sp.n. External morphology of a fixed specimen, shown 
in a picture (B) and diagram (C). 
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Figure 3. A) Spermathecae in segments 9 and 10. Note the double spermathecae in segment 10 of this 
specimen. B) Nephridial bladder of segment 7. 



A posterior gizzard is not well differentiated. There is a slight dilatation of the 
oesophagus between 14 and 16, but it lacks the muscular wall and reinforcements of a 
true gizzard. First section of the intestine is not dilated. 

Typhlosole begins around segments 20 and 21 with seven lamellae, which around 
segments 26-27 increase to nine. From there they decrease gradually in number until 
segments 80-105, where they fuse in a single lamella. The latter extends until segments 
218-230, where the typhlosole ends. 

Fraying testes and iridescent seminal funnels in 1 0 and 1 1 . Two pairs of voluminous, 
grainy seminal vesicles in 1 1 and 12. Ovaries and female funnels in 13, ovisacs in 14. 

Two pairs of spermathecae in intersegments 9/10 and 10/11 (but apparently con- 
tained in segment 9 due to septum's backward displacement), the posterior pair bigger. 
They are sessile and disc-shaped, with multiple inner chambers which open to the exte- 
rior through a common pore, in the intersegments 9/10 and 10/11. Some individuals 
show double spermathecae (each multicameral and with own pore), either in 9/10 or 
10/11 (Figure 3a). 

Anterior nephridial bladders U-shaped with very close branches and no apparent 
cecum (Figure 3b). Bladders gradually flatten towards the end of the body, taking the 
usual elongated shape. 

Distribution. Known only from its type locality. 

Habitat. The specimens were collected at 10-20 cm deep in the soil in a cleared 
holm-oak wood, at the border between a dense forest of Quercus rotundifolia and a 
dryland farm. The soil had the following characteristics: 23.03% coarse sand, 8.06% 
fine sand, 5.33% coarse silt, 60.74% fine silt, and 2.84% clay, constituting a silty loam 
soil, carbon (C): 2.40%, nitrogen (N): 0.24%, C/N: 10.18, pH: 7.98 . Mean annual 
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- Aporrectodea trapczoidcs 



Lumbricus terrestris 



- Amynthas robustus 



— Vignysa vedovinii 
Vignysa popi 

Hormogaster samnittca lirapora 

Hormogaster red/7 insularis 

Hormogaster redii 



Hormogaster redii 

■ Hormogaster redii 



- Hormogaster abbatissae 

- Hormogaster sylvestris 
Hormogaster pretiosa nigra 

- Hormogaster catalaunensis 
Hormogaster gallica 

Hormogaster arenicola 

it 1 Hormogaster riojana 

* Hormogaster ireguana 

Hormogaster eserana 
Hormogaster pretiosa var. Peralba 
Hormogaster huescana 

Hormogaster pretiosiformis 
Hormogaster oroeli 

■k I Hormogaster pretiosa arrufati 

— I Hormogaster sp. Talarn 

— Hormogaster sp. Cervera del Maestre 
i — Hormogaster castillana 
Hemigastrodrilus monicae 



*i — Hormi 
~\±r—Hc 
I * 



- Hemigastrodrilus monicae 
- Hormogaster pretiosa Villamassargia 



- Hormogaster sp. HPA 

Hormogaster najaformis 

-Xana omo deoi 

1 Hormogaster joseantonioi \ 

1/100 — EPAR1 

EBOA1 



- Dichogaster saiiens 
■ Pontodrilus litoralis 




Hormogaster elisae 
species complex 



0.2 

Figure 4. Bayesian inference of the phylogenetic tree on the concatenated sequence. Numbers above 
branches indicate posterior probability/bootstrap (of the Maximum Likelihood analysis) support values 
higher than 0.9/70 (shown as asterisks on terminal branches). Black rectangles show clades not recovered 
in both analyses (the alternative is shown with a dashed line). The cryptic species included in H. elisae are 
numbered from 1 to 5(following Novo et al. 2010). 



temperature is 12.7 °C and mean annual precipitation is 447.2 mm, as indicated by 
the nearest weather station (in Herrera de Los Navarros, Zaragoza-23 km away http:// 
www.aragon.es/DepartamentosOrganismosPublicos/Organismos/InstitutoAragon- 
esEstadistica/AreasTematicas/l4_Medio_Ambiente_Y_Energia/ci.05_Clima_Datos_ 
climatologicos.detalleDepartamento?channelSelected=ea9fa856c66de310VgnVCM2 
000002f5 5 1 bacRCRD#section 1 ) . 
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Etymology. The species is named after Jose Antonio Fernandez Fernandez, father 
of the first author Daniel Fernandez Marchan and important contributor during the 
sampling campaign in which this species was discovered. 

Molecular characters. Analyses were conducted on sequences from loci COI (10 
individuals), 16S (2 individuals), 28S (2 individuals) and H3 (2 individuals) of the 
new species, combined with similar sequences from other hormogastrid species. 

The resulting Bayesian inference of the phylogenetic tree is shown in Figure 4. Its 
topology was congruent with that of the Maximum Likelihood inferred tree, except for 
the different placement of X. omodeoi. H. joseantonioi sp.n. was recovered as a mono- 
phyletic clade, with the H. elisae species complex as a sister clade. 

Uncorrected pairwise distances for the genes COI and 16S-tRNA for H. joseanto- 
nioi and the species within the same clade (with H. elisae divided into its five cryptic 
species) are shown in Table 2. 



Discussion 

Both morphological and molecular characters of H. joseantonioi sp.n. separate it clearly 
from all known hormogastrid species, the number of typhosole lamellae and the kind 
and location of the spermathecae being particularly distinctive. Those characters, while 
failing to resolve internal relationships within Hormogastridae, have been shown to be 
suitable for species diagnosis (Rota, 1993 on typhlosole importance; Novo et al. 2012 
on spermathecae number to separate H. abbatissae from H. sylvestris). 

The species H. riojana, while distantly related according to molecular phylogeny, 
shows many similarities in morphology to H. joseantonioi (Table 3). However, H. jose- 
antonioi differ by its lower number of lamellae in its typhlosole and shorter tubercula 
pubertatis. Moreover it is longer and heavier. While the two species share a very similar 
position and shape of the spermathecae, some H. joseantonioi individuals show an ad- 
ditional spermatheca in segment 10 (on the right or left side). These cases don't seem 
to be teratologic, as the supernumerary spermathecae have their own pore in the body 
surface and contain sperm, thus being fully functional. 

Other hormogastrid species possess double or multiple spermathecae, but never of 
the multicameral, disc shaped kind. 

The geographically closest species, H. castillana (from Puerto Querol, Castellon), is 
neither morphologically nor phylogenetically closely related (Table 3). 

H. joseantonioi sp. n. appears nested on a weakly supported clade on the phyloge- 
netic tree, consisting in Hemigastrodrilus monicae, Xana omodeoi, Hormogaster pretiosa 
from Villamassargia, Hormogaster najaformis (and HPA from Omodeo, see Novo et al. 
201 1) and Hormogaster elisae. Genetic distances were high in all cases (16.47-19.08% 
for COI, 12.50-17.40% for 16S) according to the reference intervals given by Chang 
and James (201 1). Aside from H. elisae, none of them showed significant morphologi- 
cal likeness to the new species, with the very different spermathecae configurations be- 
ing noteworthy (Table 3). 
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Table 3. Comparison of the morphological characters of H. joseantonioi sp. n. and some of the phyloge- 
netically closest species (H. elisae, X. omodeoi and H. najaformis Qiu & Bouche, 1998) plus the distantly 
related H. riojana and H. castillana Qiu & Bouche, 1998. N. segments: number of segments. N. typhlo- 
sole lamellae: number of typhlosole lamellae. Body length, weight and number of segments refer to adult 
specimens. 





H. joseantonioi 


//. elisae 


X. omodeoi 


//. najaformis 


H.riojana 


//. castillana 


Colour 


Brownish 


Colourless 


Colourless 


Slightly greyish 


Dark brownish 


Brownish grey 


Clitellum 


(13)14-28 


(13)14(15)- 
26(27)28 


14-26 


13-31 


13,14,17-27,28 


1/14,15-29,1/2 30 


Tubercula 
pubertatis 


1/n 22-27 (1/n 
28) 


22(23)- 
25(26) 


23-26 


20-26 


(20)21-27 


22-28 


Length (mm) 


178-180 


92-200 


20-161 


188-230 


154 


200-325 


N. segments 


305-369 


205-300 


190-230 


395-523 


243-278 


320^29 


Weight (g) 


7.05-11.57 


1.96-9.67 


0.59^.23 


22.6-31.4 


6.57 


12.85-29.38 




9 (see text) 


9,10 
(9/10,10/11) 
Simple 
Tubular 


10,11 
(9/10,10/11) 
Simple Small, 
globular 


10,11 
(10/11,11/12) 

Multiple 
Small, globular 


9,10 




Spermathecae 
position (pores) 
and appearance 


(9/10,10/11) 
Simple(double) 
Multicameral, 

disc shaped 


(9/10,10/11) 

Simple 
Multicameral, 
disc shaped 


9,10(9/10,10/11) 
Simple Globular 


N. typhlosole 
lamellae 


9 


5 


12 


15-17 


15 


21-23 


Thickened 


6/7,7/8,8/9, 


6/7,7/8,8/9, 


(6/7),7/8,8/9, 


6/7,7/8,8/9, 


7/8,8/9,9/10, 


7/8,8/9,9/10, 


septa 


(9/10) 


(9/10) 


9/10,(10/11) 


(9/10) 


(10/11) 


(10/11) 



The H. elisae morphospecies was recovered as sister clade to H. joseantonioi sp. n. 
with high support. From a morphological point of view, most of their external char- 
acters overlap, except for a slightly longer clitellum and tubercula pubertatis, bigger 
average size and stronger pigmentation in H. joseantonioi sp. n. However, internal 
characters are very different and these species match neither in the number of lamellae 
in the typhlosole (five versus nine) nor in the structure of the spermathecae, which are 
tubular in H. elisae and disc-shaped and multicameral in H. joseantonioi. It's worth 
noting that H. elisae shares the backwardly displaced disposition of the 9/10 septum. 

Based on their phylogenetic and morphological relatedness, an origin of H. elisae 
from a common ancestor with H. joseantonioi sp. n. seems likely. This scenario is sensi- 
ble from a biogeographical point of view, as the locality of the new species is intermedi- 
ate between the ranges of H. elisae and the northeastern main hormogastrid range. A 
connection of emerged lands would have been possible from the Cretaceous-Tertiary 
boundary onwards (Andeweg 2002). 

While H. joseantonioi status as a good species and its phylogenetic relationships 
seem quite clear, generic assignment is a more problematic matter. Novo et al. (201 1) 
recovered the genus Hormogaster as paraphyletic in their molecular phylogeny, point- 
ing out the need for a deep taxonomical revision of the family Hormogastridae, cur- 
rently in preparation (author's work in progress). 

Based on its distinctive morphology and geographic range, high genetic divergence 
and consistent recovery as a well-defined clade, Novo (2010) suggested the H. elisae 
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species complex should be established as an independent genus. Due to the close phy- 
logenetic position and morphological similarity of H. joseantonioi to this clade it could 
be argued they both should be included in the same genus. 

At this stage it is more conservative to assign H. joseantonioi to the genus Hormogaster 
until the revision of the family is completed, which will allow to establish (if possible) a 
well-founded genera system on Hormogastridae. This work narrows the discontinuity 
between the North-Eastern and Central ranges of the Spanish hormogastrids. At the 
same time it highlights the importance of an intensive sampling of the area between 
Teruel and the center of the Iberian Peninsula (mainly zones of Soria and Guadalajara) 
to hopefully find new species along the hypothetical colonization route. 
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Appendix 



Supplementary material. GenBank accession numbers for all sequences used in the phylogenetic 
analysis, including outgroups. RF: sequences provided by Rosa Fernandez. 



Species 


COI 


16S-tRNAs 


28S-rRNA 


H3 


H. castillana QUE 


HQ621989 


HQ621883 


HQ62 1960.1 


HQ622028 


H. elisae 3 ANC 


EF653870 


GQ409754.1 


GQ409657.1 


HQ622001 


H. elisae 4 BOA 


GQ409661.1 


GQ409704.1 


GQ409656.1 


HQ622004 


H. elisae 1 CAB 


GQ409689.1 


GQ409729.1 


GQ409653.1 


HQ622007 


H. elisae 1 FRE 


GQ409698.1 


GQ409723.1 


GQ409653.1 


HQ622009 


H. elisae 1 JAR 


GQ409665.1 


GQ409745.1 


GQ409653.1 


HQ622013 


H. elisae 1 LOZ 


EF653888 


GQ409725.1 


GQ409653.1 


HQ622016 


H. elisae 1 MOL 


EF653875 


GQ409732.1 


GQ409653.1 


HQ622019 


H. elisae 1 NAV 


GQ409683.1 


GQ409730.1 


GQ409653.1 


HQ622021 


H. elisae 5 PAR 


EF653898 


GQ409709.1 


GQ409655.1 


HQ622024 


H. elisae 1 RED 


EF653881 


GQ409741.1 


GQ409653.1 


HQ622029 


H. elisae 4 SEV 


EF653905 


GQ409707.1 


GQ409656.1 


HQ622031 


H. elisae 2 SIG 


EF653893 


GQ409710.1 


GQ409654.1 


HQ622033 


H. elisae 2 SOT 


GQ409700.1 


GQ409716.1 


GQ409654.1 


HQ622034 


H. elisae 1 TRE 


GQ409678.1 


GQ409737.1 


GQ409653.1 


HQ622038 


H. elisae 1 UCE 


GQ409692.1 


GQ409720.1 


GQ409653.1 


HQ622039 


H. elisae 1 VEN 


GQ409671.1 


GQ409750.1 


GQ409653.1 


HQ622041 


H. pretiosa arrufati 


HQ621995 


HQ621889 


HQ62 1966.1 


HQ622040 


H. pretiosa var. PRB 


HQ621987 


HQ621881 


HQ621958.1 


HQ622026 


H. pretiosa Villamassargia 


HQ621998 


HQ621893 


HQ62 1969.1 


HQ622045 


H. pretiosiformis oroeli 


HQ621984 


HQ621877 


HQ621955.1 


HQ622022 


H. redii redii 


HQ621978 


HQ621871 


HQ62 1949.1 


HQ622012 


H. redii redii 


HQ621971 


HQ621863 


HQ62 1942.1 


HQ622000 


H. redii redii 


HQ621976 


HQ621869 


HQ62 1947.1 


HQ622010 


H. redii insularis 


HQ621996 


HQ621890 


HQ62 1967.1 


HQ622042 


H. samnitica lirapora 


HQ621993 


HQ621887 


HQ62 1964.1 


HQ622036 


Hemigastrodrilus monicae 


HQ621979 


HQ621872 


HQ62 1950.1 


HQ622014 


Hemigastrodrilus monicae 


HQ621982 


HQ621875 


HQ621953.1 


HQ622018 


Hormogaster abbatissae 


HQ621990 


HQ621884 


HQ621961.1 


HQ622030 


Hormogaster arenicola 


HQ621972 


HQ621865 


HQ62 1943.1 


HQ622003 


Hormogaster catalaunensis 


HQ621973 


HQ621866 


HQ62 1944.1 


HQ622005 


Hormogaster eserana 


HQ621977 


HQ621870 


HQ62 1948.1 


HQ622011 


Hormogaster gallica 


HQ621974 


HQ621867 


HQ621945.1 


HQ622006 


Hormogaster huescana 


HQ621980 


HQ621873 


HQ621951.1 


HQ622015 


Hormogaster ireguana 


HQ621994 


HQ621888 


HQ621965.1 


HQ622037 


Hormogaster najaformis 


HQ621985 


HQ621878 


HQ621956.1 


HQ622023 


Hormogaster nigra 


HQ621988 


HQ621882 


HQ621959.1 


HQ622027 


Hormogaster pretiosiformis 


HQ621983 


HQ621876 


HQ621954.1 


HQ622020 


Hormogaster riojana 


HQ621970 


HQ621862 


HQ621941.1 


HQ621999 


Hormogaster sp. CER 


HQ621975 


HQ621868 


HQ62 1946.1 


HQ622008 


Hormogaster sp. HPA 




HQ621892 




HQ622044 


Hormogaster sp. TAL 


HQ621992 


HQ621886 


HQ621963.1 


HQ622035 
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Species 


COI 


16S-tRNAs 


28S-rRNA 


H3 


Hormogaster sylvestris 


HQ621981 


HQ621874 


HQ621952.1 


HQ622017 


Vignysa popi 


HQ621991 


HQ621885 


HQ62 1962.1 


HQ622032 


Vignysa vedovinii 


HQ621986 


HQ621880 


HQ621957.1 


HQ622025 


Xana omodeoi 


HQ621997 


HQ621891 


HQ62 1968.1 


HQ622043 


Amynthas robustus 


EF077569.1 


EF490524.1 


EF490529.1 




Dkhogasters aliens 




AF406573.1 


AY101560.1 




Pontodrilus litomlis 




AY340473.1 






Lumbricus terrestris 


HQ691222 


U24570 


HQ691218 


HQ691227 


Aporrectodea trapezoides 


RF 


HQ621864 


RF 


HQ622002 



